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ABSTRACT 

Aims. We present an asteroseismological study on PG 0122+200, the coolest known pulsating PG1159 (GW Vir) star. 
Our results are based on an augmented set of the full PG1159 evolutionary models recently presented by Miller Bertolami 
& Althaus (2006). 

Methods. We perform extensive computations of adiabatic g-mode pulsation periods on PG1159 evolutionary models 
with stellar masses ranging from 0.530 to O.741M0 . These models take into account the complete evolution of progenitor 
stars, through the thermally pulsing AGB phase and born-again episode. We constrain the stellar mass of PG 0122+200 
by comparing the observed period spacing with the asymptotic period spacing and with the average of the computed 
period spacings. We also employed the individual observed periods to find a representative seismological model for 
PG 0122+200. 

Results. We derive a stellar mass of 0.626Mq from a comparison between the observed period spacing and the computed 
asymptotic period spacing, and a stellar mass of O.567M0 by comparing the observed period spacing with the average 
of the computed period spacing. We also find, on the basis of a period-fit procedure, an asteroseismological model 
representativ e of PG 0122+200 which is able to reproduce the observed period pattern with an average of the period 
differences of (5IL = 0.88 s and a root-mean-square residual of a sn . = 1.27 s. The model has an effective temperature 
T e s = 81500 K, a stellar mass M* = O.556M0, a surface gravity logy = 7.65, a stellar luminosity and radius of 
log(L*/Z/0) = 1.14 and \og(R t / Rq) = —1.73, respectively, and a He-rich envelope thickness of M cnv = 1.9 x 1O _2 M0. 
We derive a seismic distance d ~ 614 pc and a parallax 7r ~ 1.6 mas. The results of the period- fit analysis carried 
out in this work suggest that the asteroseismological mass of PG 0122+200 could be ~ 6 — 20% lower than thought 
hitherto and in closer agreement (to within ~ 5%) with the spectroscopic mass. This result suggests that a reasonable 
consistency between the stellar mass values obtained from spectroscopy and asteroseismology can be expected when 
detailed PG1159 evolutionary models are considered. 

Key words, stars: evolution — stars: interiors — stars: oscillations — stars: variables: other (GW Virginis) — white 
dwarfs 



1. Introduction 

PG 0122+200 (BB Psc or WD 0122+200) is the coolest 
known pulsating PG1159 star belonging to the GW Vir 
class of variables. GW Vir stars are very hot hydrogen- 
deficient post- Asymptotic Giant Branch (AGB) stars with 
surface layers rich in helium, carbon and oxygen (Werner & 
Herwig 2006). They exhibit multiperiodic luminosity varia- 
tions with periods in the range 5 — 50 minutes, attributable 
to nonradial pulsation g-modes. PG1159 stars are thought 
to be the evolutionary link between Wolf-Rayet type cen- 
tral stars of planetary nebulae and most of the hydrogen- 
deficient white dwarfs. It is generally accepted that these 
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stars have their origin in a born-again episode induced by 
a post-AGB helium thermal pulse (see Iben et al. 1983, 
Herwig et al. 1999, Lawlor & MacDonald 2003, Althaus et 
al. 2005, Miller Bertolami et al. 2006 for recent references). 

PG 0122+200 is characterized by T e $ = 80 000±4 000 K 
and \ogg = 7.5 ± 0.5 (Dreizler & Heber 1998). At this 
effective temperature, PG 0122+200 currently defines the 
locus of the low-luminosity red edge of the GW Vir insta- 
bility strip. The photometric variations of this star were 
discovered by Bond & Grauer (1987). Besides the intrin- 
sic interest in probing its interior, pulsation studies of 
PG 0122+200 offer a unique opportunity to study neutrino 
physics. Indeed, at the evolutionary stage characterizing 
PG 0122+200, neutrino emission constitutes a main energy 
sink 1 (O'Brien et al. 1998). 



At variance with the solar neutrino emission, which is a by- 
product of nuclear fusion, the neutrino flux of pre-white dwarf 
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The determination of the stellar mass of PG 0122+200 
has been the subject of numerous investigations. The stel- 
lar mass of pulsating pre-white dwarfs can be constrained, 
in principle, from asteroseismology — the asteroseismologi- 
cal mass — either through the observed period spacing (see, 
for instance, Kawaler & Bradley 1994; Corsico & Althaus 
2006) or by means of the individual observed periods (see, 
e.g., Kawaler & Bradley 1994, Corsico & Althaus 2006, 
Corsico et al. 2007). The early study of O'Brien et al. 
(1996) predicts a stellar mass of about 0.66 — O.72M for 
PG 0122+200 corresponding to an observed mean period 
spacing of 21.2 s. Vauclair et al. (1995), on the other hand, 
suggest an even higher stellar mass, based on a observed 
mean period spacing of PG 0122+200 of ~ 16 s. Later, 
O'Brien et al. (1998) find strong evidence for a I = 1 mean 
period spacing of 21 s, although a value of ~ 16 s can- 
not be conclusively ruled out. These values of the period 
spacing imply a stellar mass of ~ 0.69M Q and ~ l.OM , 
respectively, based on the PG 1159 models then available. 
By means of a period- fit procedure based on PG1159 evo- 
lutionary models with several masses derived from the full 
sequence of O.589M of Althaus et al. (2005), Corsico & 
Althaus (2006) obtain a stellar mass of M* = O.64M 
for PG 0122+200. Recently, Fu et al. (2007) (hereinafter 
FUEA07) have presented new multisite photometric obser- 
vations of PG 0122+200 obtained in 2001 and 2002. By col- 
lecting the new data together with previous observations, 
these authors have succeeded in detecting a total of 23 fre- 
quencies corresponding to modes with £ = 1 and derived 
unambiguously a mean period spacing of 22.9 s. On the 
basis of the models of Kawaler & Bradley (1994), these au- 
thors inferred a stellar mass of 0.59 ± O.O2Af . 

The stellar mass of PG1159 stars can also be esti- 
mated through the comparison of the spectroscopic values 
of T e ff and g with evolutionary tracks — the spectroscopic 
mass. On the basis of the evolutionary tracks of O'Brien 
& Kawaler (2000), Dreizler & Heber (1998) derived a stel- 
lar mass of 0.53 ± 0.1M© for PG 0122+200. On the other 
hand, Werner & Herwig (2006) determined M* = O.58M 
from a comparison with the H-rich evolutionary models of 
Schonberner (1983). The most recent determination is that 
of Miller Bertolami & Althaus (2006), who derived a stellar 
mass of 0.53Mq on the basis of PG1159 evolutionary mod- 
els that take fully into account the evolutionary history and 
the surface composition of the progenitor stars. 

The discrepancy between the asteroseismological mass 
derived by FUEA07 (0.57 < M*/M < 0.61) and the most 
recent spectroscopic determination (O.53M ) has prompted 
us to undertake the present asteroseismological investiga- 
tion for PG 0122+200, taking full advantage of the new 
generation of PG1159 evolutionary models recently devel- 
oped by Miller Bertolami & Althaus (2006). These authors 
have followed in detail all of the evolutionary phases prior to 
the formation of PG1159 stars with different stellar masses, 
particularly the born-again stage. In addition to the issue of 
the stellar mass, the employment of such detailed PG1159 
models allows us to address the question of the He-rich 
envelope mass (q v = M cnv /M Q ) of PG 0122+200, which 
FUEA07 constrain to be in the range —6 <^ iogq y <J —5.3. 



stars such as PG 0122+200 is the result of different scattering 
processes, being the plasmaneutrino, Bremsstrahlung neutrino 
and photoneutrino emission the most relevant ones (see O'Brien 
& Kawaler 2000). 
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Fig. 1. Our PG1159 full evolutionary tracks in the logT c ff — 
log g plane, labelled with the corresponding stellar mass value 
in solar units. The black circle is the location of PG 0122+200 
according to spectroscopy (T c g = 80 ± 4 kK and log g — 7.5 ± 
0.5). Note the large error box (dashed), in particular for logg. 
The square (red) is the location of the star as predicted by our 
asteroseismological analysis (see §4). The blue (hot) boundary of 
the theoretical dipole (£ — 1) instability domain — according to 
Corsico et al. (2006) — is also shown [Color figure only available 
in the electronic version of the article.] 



Finally, a precise knowledge of the mass of PG 0122+200 
is a crucial aspect concerning the role played by neutrinos 
in that star. 

The paper is organized as follows: in the next Section we 
briefly describe our PG1159 evolutionary models. In Sect. 3 
we derive the stellar mass of PG 0122+200 by means of the 
observed period spacing. In Sect. 4 we derive structural pa- 
rameters of this star by employing the individual observed 
periods. In this section we derive an asteroseismological 
model representative of PG 0122+200 (§ 4.1) and discuss 
its main structural and pulsational characteristics (§ 4.2), 
its helium envelope thickness (§ 4.3), its mode-trapping 
properties (§ 4.4) and the asteroseismological distance and 
parallax (§ 4.5). Finally, in Sect. 5 we summarize our main 
results and make some concluding remarks. 



2. Evolutionary models and numerical tools 

The pulsation analysis presented in this work relies on a 
new generation of stellar models that take into account 
the complete evolution of PG1159 progenitor stars. These 
models have been recently employed by our group for 
a pulsation stability analysis of the GW Vir stars and 
for an asteroseismological study of the hot PG1159 star 
RX J2117.1+3412 (Corsico et al. 2006 and Corsico et al. 
2007, respectively). The stages for the formation and evo- 
lution of PG1159 stars were computed with the LPCODE 
evolutionary code, which is described at length in Althaus 
et al. (2005). The neutrino production rates adopted in our 
computations are those of Itoh et al. (1989, 1992). 



Corsico et al.: Asteroseismological constraints on PG 0122+200 



3 




Fig. 2. The dipole (£ = 1) asymptotic period spacing (All?) 
for different stellar masses in terms of the effective temperature. 
Numbers along each curve denote the stellar masses (in solar 
units). The plot also shows the location of PG 0122+200 (T cff = 
80 + 4 kK and All = 22.9 s) and the remainder high-gravity, 
low-luminosity GW Vir stars (PG 1159-035, PG 2131+066, 
and PG 1707+427) with the period spacing and T B g data taken 
from Kawaler et al. (2004). The mass of PG 0122+200 as de- 
rived by comparing An? with All is M, = 0.625T.Hl6 M © 
[Color figure only available in the electronic version of the arti- 
cle.]. 



Specifically, the background of stellar models was ex- 
tracted from the evolutionary calculations recently pre- 
sented by Althaus et al. (2005), Miller Bertolami & Althaus 
(2006), and Corsico et al. (2006), who computed the com- 
plete evolution of model star sequences with initial masses 
on the ZAMS in the range 1 — 3.75M Q . We refer the reader 
to those works for details. Suffice it to mention that all of 
the post-AGB evolutionary sequences have been followed 
through the very late thermal pulse (VLTP) and the re- 
sulting born-again episode that give rise to the H-dcficient, 
helium-, carbon- and oxygen- rich composition character- 
istic of PG1159 stars. The masses of the resulting rem- 
nants are 0.530, 0.542, 0.556, 0.565, 0.589, 0.609, 0.664, and 
0.741M Q . The new sequence with M — O.556M , coming 
from a progenitor star with M* = 1.8M at the ZAMS, has 
been computed specifically for the present asteroseismolog- 
ical study. The evolutionary tracks in the logT c ff — log g 
plane for the PG1159 regime are displayed in Fig. 1. 

It is worth mentioning that the use of these evolution- 
ary tracks constitutes a major improvement with respect 
to previous asteroseismological studies. As mentioned, our 
PG1159 evolutionary sequences are derived from the com- 
plete born-again evolution of progenitor stars and a careful 
treatment of the mixing and extramixing processes dur- 
ing the core helium burning, fundamental aspects when 
attempts are made at constructing stellar models appro- 
priate for PG1159 stars. In particular, these evolutionary 
calculations reproduce: (1) the spread in surface chemical 
composition observed in PG1159 stars, (2) the short born- 
again times of V4334 Sgr (see Miller Bertolami et al. 2006 
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Fig. 3. Same as Fig. 2, but for the average of the computed 
period spacings (All*.). The mass of PG 0122+200 as derived 
by comparing AII fe with All is M, = 0.567t° °oilM Q [Color 
figure only available in the electronic version of the article]. 

and Miller Bertolami & Althaus 2007a), (3) the location 
of the GW Vir instability strip in the log T c g — log g plane 
(Corsico et al. 2006), and (4) the expansion age of the plan- 
etary nebula of RX J2117. 1+3412 (see the paper by Corsico 
et al. 2007 and its associated erratum). We believe that the 
employment of these new evolutionary computations render 
reliability to our pulsational inferences for PG 0122+200. 

We computed £ — 1 <?-mode adiabatic pulsation pe- 
riods and asymptotic period spacings with the same nu- 
merical code we employed in our previous works (see, e.g., 
Corsico & Althaus 2006 for details). We analyzed about 
3000 PG1159 models covering a wide range of effective 
temperatures and luminosities (5.4 ;> log(T e g) J> 4.8 and 
<; log(i*/i>0) <; 4.2, respectivley) and a range of stellar 
masses (0.530 < M*/M < 0.741). 

3. Mass determination from the observed period 
spacing 

Here, we constrain the stellar mass of PG 0122+200 by 
comparing the asymptotic period spacing, An|, and the 
average of the computed period spacings, An^ (k being 
the radial order), with the observed period spacing, ALT 02 
These methods take full advantage of the fact that the pe- 
riod spacing of PG1159 pulsators depends primarily on the 
stellar mass, and weakly on the luminosity and the He-rich 
envelope mass fraction (Kawaler & Bradley 1994; Corsico 
& Althaus 2006). Note that in these approachs we make 
use of the spectroscopic constraint that the effective tem- 
perature of PG 0122+200 is of 80 kK (Dreizler & Heber 
1998). 

The asymptotic period spacing and the average of the 
computed period spacings for £ = 1 modes as a func- 

2 Note that most asteroseismological studies rely on the 
asymptotic period spacing to infer the stellar mass of GW Vir 
pulsators. 
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tion of the effective temperature are displayed in Figs. 2 
and 3, respectively, for different stellar masses. Also shown 
in these diagrams is the location of PG 0122+200, with 
An° = 22.9 s (FUEA07). Here, AH"f = U /y/e(£+l), 
where n = 2n 2 [J^ (iV/r)dr] _1 , being N the Brunt- Vaisala 

frequency (Tassoul et al. 1990). The quantity Allfe, on the 
other hand, is assessed by averaging the computed forward 
period spacings (An* = H-k+i —11*) in the range of the 
observed periods in PG 0122+200 (330-620 s; see Table 1). 
From a comparison between All and AII| we obtain 

a stellar mass of M* = 0.625to ow M e- The quoted uncer- 
tainties in the value of Af* come from the errors in the spec- 
troscopic determination of the effective temperature. In the 
same way, we get A/* = O.5671q;o? 3 M0 if we compare All 
and Allfc. The higher value of Af* (about 10% larger) as de- 
rived from AII^ is due to that usually the asymptotic period 
spacing is larger than the average of the computed period 
spacings (see Corsico & Althaus 2006), in particular for the 
short periods like those exhibited by PG 0122+200, i.e. for 
which the full asymptotic regime of the modes (fc 3> 1) 
has not been attained 3 . It is important to note that the 
first method to derive the stellar mass is somewhat less 
realistic than the second one, because the asymptotic pre- 
dictions are, in principle, only valid for chemically homoge- 
neous stellar models, while our PG1159 models are indeed 
chemically stratified. 

Finally, we note that our inferred stellar mass values 
of Af* w 0.57Af Q and in particular Af* ps 0.63M q are in 
conflict with the value Af* = 0.53Af Q as derived from spec- 
troscopy coupled to evolutionary tracks (Dreizler & Heber 
1998; Miller Bertolami & Althaus 2006). 



4. Constraints from the individual observed periods 

4.1. The search for the best- fit model 

In this approach we seek a pulsation model that best 
matches the individual pulsation periods of PG 0122+200. 
We assume that all of the observed periods correspond to 
£ = 1 modes (see FUEA07). The goodness of the match 
between the theoretical pulsation periods (life) and the ob- 
served individual periods (IIP) is measured by means of a 
quality function defined as x 2 (A/*, T e g) — Y^i=i min[(IiP — 
IIfc) 2 ]/n, where n (= 9) is the number of observed peri- 
ods (first column in Table 1). The PG 1159 model that 
shows the lowest value of x 2 will be adopted as the "best- 
fit model" . This approach has also been used by Corsico & 
Althaus (2006) and Corsico et al. (2007). 

We evaluate the function x 2 (M„, T c ff ) for stellar 
masses of 0.530,0.542,0.556,0.565,0.589,0.609,0.664, and 
0.741M Q . For the effective temperature we employed a 
much more finer grid (AT e g = 10 — 30 K). The quan- 
tity (x 2 )~* in terms of the effective temperature for dif- 
ferent stellar masses is shown in the mosaic of Fig. 4 
together with the spectroscopic effective temperature of 
PG 0122+200. We find one strong maximum for the model 
with Af* = O.556Af and T cB « 81.5 kK (panel c). Such a 



3 At variance with this, for the longer periods exhibited by 
RX J2117.1+3412 (with 30 < fc < 53) the asymptotic conditions 
are more nearly reached and, as a result, the stellar mass derived 
from the asymptotic period spacing is very close to that derived 
from de average of the computed period spacings (see Corsico 
et al. 2007). 
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Fig. 4. The inverse of the quality function of the period fit in 
terms of the effective temperature for the PG1159 sequences 
with different stellar masses indicated (in solar mass) in each 
panel. The grey vertical strip corresponds to the spectroscopic 
effective temperature of PG 0122+200 and its uncertainties 
(T cff = 80 000 ±4 000 K). Note the strong maximum in (x 2 ) -1 
for M* = O.556M at T cff » 81 500 K. This corresponds to our 
"best-fit" model (see text for details) [Color figure only available 
in the electronic version of the article]. 



pronounced maximum in the inverse of x 2 implies an excel- 
lent agreement between the theoretical and observed peri- 
ods. Another maximum, albeit somewhat less pronounced, 
is encountered for the model with Af* = 0.542Af Q at 
Toff ~ 87.7 kK and constitutes another acceptable as- 
teroseismological solution, in particular because its stellar 
mass is closer to the spectroscopic mass of PG 0122+200 
(0.53Af Q ). However, because the agreement between ob- 
served and theoretical periods for this model is somewhat 
poorer than for the one with M* = O.556Af , we choose 
to adopt this last model as the best-fit asteroseismologi- 
cal model. Note that our best-fit model has an effective 
temperature very close to that suggested by spectroscopy, 
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Table 1. Observed m = periods (II?) for PG 0122+200 
(taken from FUEA07), theoretical i = l,m = periods (life), 
period differences (SUi = II? — II fc ), radial orders (fc), linear 
growth rates (r}k), and rates of period change (II fe ) for the best- 
fit model. 



n? 




SUi 


k 




n fc 


N 


[b] 


M 




[io- 6 ] 


[IO" 12 s/s] 


336.68 


334.12 


2.56 


12 


1.12 


1.22 




354.85 




13 


1.74 


1.70 


380.10 


380.70 


-0.60 


14 


3.54 


1.84 


400.99 


400.35 


0.64 


15 


6.13 


1.59 




425.25 




16 


8.53 


2.24 


449.48 


448.16 


1.32 


17 


15.90 


1.71 


468.69 


469.32 


-0.63 


18 


16.30 


2.46 


494.92 


494.76 


0.16 


19 


28.64 


2.19 


517.96 


516.65 


1.31 


20 


36.50 


2.10 




539.65 




21 


32.24 


3.11 


564.28 


563.98 


0.30 


22 


61.06 


2.00 




585.80 




23 


48.17 


3.04 


611.15 


610.73 


0.42 


24 


55.99 


3.26 



well inside the error bar. A detailed comparison of the 
observed m = periods in PG 0122+200 with the theo- 
retical periods of the best-fit model is provided in Table 
1. The high quality of our period fit is quantitatively re- 
flected by the average of the absolute period differences, 
3Hi = (ELi I^D/n, where SU, = IT? - U k , an d by the 
root-mean-square residual, a sn . = \J |<5IL| 2 )/n. We ob- 
tain bW-i = 0.88 s and a sn = 1.27 s. The quality of our 
fit for PG 0122+200 is much better than that achieved by 
Corsico et al. (2007) for RX J2117.1+3412 (Ml = 1.08 
s) and those obtained by Kawaler & Bradley (1994) and 
Corsico & Althaus (2006) (W~ = 1.19 s and Jn~ = 1.79 
s, respectively) for PG 1159—035. Note that we are able to 
get a PG1159 model that nicely reproduces the period spec- 
trum observed in PG 0122+200 without artificially tuning 
the value of structural parameters such as the thickness of 
the outer envelope, the surface chemical abundances, or the 
shape of the core chemical profile which, instead, are kept 
fixed at the values predicted by the evolutionary computa- 
tions. 

Table 1 also shows the linear growth rates (%) of the 
fitted pulsation modes (fifth column) for our best-fit model, 
computed with the nonadiabatic pulsation code described 
in Corsico et al. (2006). We found that all of the fitted 
modes have positive values of rjk, implying pulsational in- 
stability, although our stability analysis predicts a band of 
unstable mode-periods (230 <; life ^ 730 s) somewhat wider 
than the interval of periods detected in PG 0122+200. 

The last column in Table 1 shows the rate of period 
change of the fitted pulsation modes. Our calculations pre- 
dict all of the pulsation periods to increase with time (life > 
0), in accordance with the decrease of the Brunt- Vaisala, fre- 
quency in the core of the model induced by cooling. Note 
that at the effective temperature of PG 0122+200, cooling 
has the largest effect on life , while gravitational contraction, 
which should result in a decrease of periods with time, be- 
comes negligible and no longer affects the pulsation periods, 
except for the case of modes trapped in the envelope (see 
§4.4). Until now, the only secure measurement of II in pre- 
white dwarf stars is that of PG 1159—035, the prototype of 
the class, for which Costa et al. (1999) obtained a positive 



value of Id = (+1.307 ± 0.003) x IO" 10 s/s for the 516 s 
period. Note that our theoretical (positive) life values for 
the best-fit model (1.22-3.26 x 10~ 12 s/s) are two orders of 
magnitude lower. For the case of PG 0122+200, a determi- 
nation of any II has not been assessed up to now, although 
work in this direction is in progress (see Fu et al. 2002). 

4.2. Characteristics of the best-fit model 

The main features of our best-fit model are summarized 
in Table 2, where we also include the parameters of 
PG 0122+200 from other published studies 4 . Note that the 
effective temperature of our best-fit model is virtually the 
same as the spectroscopic value. Thus, the location of the 
star in the logT e ff — logg plane is vertically shifted to higher 
gravities according to our predictions (see Fig. 1). 

Our best-fit model has a stellar mass of M* = 0.556M Q , 
somewhat smaller than the value derived from the aver- 
age of the computed period spacing, M* ~ 0.57M Q , and 
substantially lower than that inferred from the asymp- 
totic period spacing, M* <~ 0.63M Q (see §3). On the other 
hand, FUEA07 have inferred a value of the stellar mass of 
PG 0122+200 by using an interpolation formula to the pe- 
riod spacing derived by Kawaler & Bradley (1994) on the 
basis of a large grid of artificial PG1159 models in the lumi- 
nosity range 1.6 <J log(L*/L Q ) <j 3.0. These authors obtain 
a rather high value of 0.59 ± 0.02M Q , in line with the trend 
of early determinations (O'Brien et al. 1998) and also in 
good agreement with our values derived from the period 
spacing, but in clear conflict with the mass of our best-fit 
model. 

On the other hand, the M* value of our best-fit model 
is somewhat higher than the spectroscopic mass of 0.53M Q 
derived by Miller Bertolami & Althaus (2006) (see also 
Drcizler & Heber 1998) for PG 0122+200. Note that a dis- 
crepancy between the asteroseismological and the spectro- 
scopic values of M„ is generally encountered among PG1159 
pulsators (see Corsico et al. 2006, 2007). Until now, the 
asteroseismological mass of PG 0122+200 has been about 
10 - 30% larger (AM, w 0.06 - 0.17M Q ) than the spec- 
troscopic mass. In light of the best-fit model derived in 
this paper, this discrepancy is notably reduced to less than 
about 5% (AM, « 0.026M Q ). 

FUEA07 infer the stellar luminosity of PG 0122+200 by 
using the formula of Kawaler & Bradley (1994) mentioned 
above. They obtain log(L»/L Q ) = 1.3±0.5, larger than the 
luminosity of our best-fit model, log(L»/L ) = 1.14+q'q 4 , 
and with an accuracy a factor 20 worse. The large uncer- 
tainty in the luminosity quoted by FUEA07 is due to the 
large uncertainty in the spectroscopically determined logg, 
a quantity used by these authors to derive the luminosity. 

4.3. Helium-rich envelope thickness 

An important parameter to be discussed separately is the 
thickness of the outer envelope (M cnv ) of PG 0122+200. 
We define M env as the mass above the chemical discon- 
tinuity between the He-rich envelope and the C/O core. 

4 Errors in T e g and log(L*/L ) are estimated from the width 
of the maximum in the function \ 2 vs - ^eff and log(L*/LQ), 
respectively; the error in the stellar mass comes from the grid 
resolution in M*. Errors in the remainder quantities are derived 
from these values. 
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Table 2. The main characteristics of PG 0122+200. The second column corresponds to spectroscopic results, whereas the third 
and fourth columns present results from the pulsation study of FUEA07 and from the asteroseismological model of this work, 
respectively. 



Quantity 


Spectroscopy 




Asteroseismology 

I ± Ills WOIKI 


Jeff |Krv] 


°.n 4- a ( a ) 

oU in 4 




81.04_ 1 4 


M* [Mq\ 


n 4- n i (t>) 
U.Oo ± U.l v y 


u.oy ± U.Uz 


0.556_ 014 


log gr [cm/s 2 ] 


7.5±0.5 (a) 


— 




log(L*/L ) 


1 <^+0.2(**) 

^-0.3 


1.3 + 0.5 


1 14+0-02 


log{R,/R e ) 


-i-eslg; 1 ^"' 


-1.65 + 0.25 


i yq - I _ 0.025 

— i.(o_ 01 


M cnv [M Q ] 




(6 - 30) x 10~ 7 


0.019 + 0.006 


C/He, 0/He w 


0.9,0.4 (a) 




1.14,0.71 


BC [mag] 


_c: O1+0.23 
°- oi -0.21 




-5.89t°;°* 


My [mag] 


7 rr+0.74 




7 7Q+0.03 
' • ' y -0.10 


Mboi [mag] 


1.74 




1 9+0" 

1 - y -0.14 


A v [mag] 


0.19 




0.071 


d [pc] 


682 


700+ 1000 

' uu -400 


61411 


7T [mas] 


1.47 




1.6 + 0.1 



Note: (*) Abundances by mass, (**) Interpolated from the tracks by assuming spectroscopic (T c ff , log fit) = (80kK, 7.5). 
References: (a) Dreizler & Heber (1998); (b) Miller Bertolami & Althaus (2006). 



Our best-fit model has M env = 0.019M Q . On the other 
hand, FUEA07 suggest a value of M cnv in the range 
(6 — 30) x 10~ 7 M©, about 5 orders of magnitude smaller. 
In this section, we try to answer the question: could a strik- 
ingly low value of M cnv like that suggested by FUEA07 be 
explained by mass loss during the PG1159 phase? To this 
end, we performed additional PG1159 evolutionary calcu- 
lations to explore the amount of stellar mass that could be 
eroded by winds. Specifically, we have performed new evo- 
lutionary simulations for the sequence of the best-fit model 
(M* = 0.556 M©) starting from the second departure (post- 
VLTP) of the AGB until the PG1159 stage is reached, with 
different mass loss rate prescriptions. Specifically, we have 
adopted two different mass loss rates (Mi, M 2 ) appropri- 
ate for radiatively driven winds. Namely, the one given by 
Blocker (1995), which is based on Pauldrach et al. (1988), 
results 

/ L \ 186 

Mi = 1.29 x 10- 15 [M /yr], (1) 

and the one adopted by Lawlor & MacDonald (2006), 
which is based on a modified version of the treatment of 
Abbott (1982), 

M 2 = 1.2X10- (0' 2 [M e/yr| .( 2) 

In the last expression M c g = (1 — r)M* with V defined as 
in Castor et al. (1975). The metallicity was set to Z = Z@, 
because at high metallicities iron lines are expected to be 
dominant for radiative driven winds (Vink et al. 2001) and 
iron abundance is expected to remain unchanged during the 
whole evolution 5 . Roughly, Mi is about one order of mag- 
nitude lower than Mi in the present simulations. The total 

5 However, how important C,N,0 lines can be at the extremely 
high abundances of PG1159 stars is not known. In any case, we 
think that the inclusion of a simulation with M3 = IOM2 — that 
would correspond to the inclusion of a value Z/Zq — 100, or 
"Z = 2" in Eq. (2); see the text — really sets an upper limit for 
possible mass loss rates during the evolution. 



amount of mass lost by these sequences when they reach 
the location of PG 0122+200 is 7 x 1O~ 5 M for Mi and 
4.4 x 10~ 4 M Q for M2, which are both negligible as com- 
pared with the mass of the envelope of the best-fit model. 
For completeness we have considered a more extreme case 
by adopting a mass loss rate of M3 = 10 M2. In this case 
the mass loss rate at the WR-CSPN stage (L* - 10 000 L Q 
and T cff < 100 000 K) is of the order of several 10" 6 M /yr, 
and the rate at the evolutionary "knee" in the HR dia- 
gram during the PG1159 stage is of about 10~ 7 M Q /yr. 
These values are consistent with the largest rates observed 
at both PG1159 and WR-CSPN stages (Koesterke et al. 
1998, Koesterke 2001) and, consequently, are probably an 
overestimation of the effect in view of the low mass of our 
best-fit model. Even in this case, the mass eroded by winds 
amounts to only 3.4 x 10 -3 M© which is about one order 
of magnitude lower than the initial mass of the envelope 6 . 
Thus, it seems that envelopes as thin as those proposed 
by FUEA07 could be ruled out in the context of single 
star stellar evolution. More importantly, the reduction in 
the mass of the He-rich envelope from a canonical value 
of ~ 1O~ 2 M to a value of ~ 10~ 7 M Q would require an 
extreme fine tuning (of five orders of magnitude) in the 
mass-loss rate to avoid the complete removal of the whole 
envelope. In the absence of a mechanism that justifies this 
fine tuning, such extremely thin envelopes should be taken 
with some caution. 

4.4. Mode trapping 

In this section we shall try to disentangle the possi- 
ble mode-trapping signatures that could be hidden in 
the observed period spectrum of PG 0122+200. Following 
FUEA07, we consider the residuals (Rn) of the period dis- 
tribution relative to the mean period spacing 7 . For the 

6 It is interesting to note that, even for this extreme case, the 
period-fit does not deteriorate significantly as compared with 
the case of the best-fit model (see Sect. 4.1). 

7 Residuals relative to the mean period spacing are more ap- 
propriate than the forward period spacing (AITfc = rU+i — life) 
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Fig. 5. Panel a: Distribution of the residuals Rn relative to 
the mean period spacing for the case of the observed periods 
(red) and for the case of the calculated periods (black) of the 
best-fit model. Panel b: the distribution of the kinetic energy. 
Panel c: the values of the relative rates of period change. Panel 
d: the values of the linear nonadiabatic stability coefficients 
77 = — 9(a)/K(cr) (being a the complex eigenfrequency) . The 
numbers correspond to the radial order k of the modes trapped 
in the envelope. See the text for details [Color figure only avail- 
able in the electronic version of the article! . 



case of PG 0122+200, a linear least-square fit to the ob- 
served periods gives a mean period spacing of 22.97 s, 
while for the best-fit model we obtain a mean theoreti- 
cal period spacing of 23.05 s. In panel a of Fig. 5 we plot 
the i?n-distribution for the observed periods (red) and for 
the case of the calculated periods (black) of the best fit 
model. The calculated distribution is in very good agree- 
ment with the observed one, in particular for the modes 
with k = 15, 17, 19, 20, 22, 24. In addition, the global struc- 
ture of maxima and minima seen in the observed distribu- 
tion is nearly duplicated by the computed one. 

Mode trapping in PG1159 stars has been discussed at 
length by Kawaler & Bradley (1994) and Corsico & Althaus 
(2006); we refer the reader to those works for details. Here, 
we shall try to answer the question: which modes could be 
trapped in the outer envelope of PG 0122+200? At first 
glance, they would be the modes showing a minimum in 



the i?n — n diagram. As we shall see below, this criterion 
can lead to erroneous conclusions. A more secure way to 
find which modes are trapped in the outer envelope is to 
examine their pulsation kinetic energy (i?kin)- I n panel b 
of Fig. 5 we show the kinetic energy distribution for our 
best- fit model 8 . Since modes that oscillate mainly in the 
outer envelope have lower Skin values, one can easily iden- 
tify trapped modes as those having local minima in the 
kinetic energy distribution. As can be seen from the figure, 
they are the modes with k = 12, 15, 17, 20 and 22. Note 
that in some cases a minimum in i?n does coincide with a 
minimum in log_Ekin (for instance for k = 15) and in other 
cases does not (for instance for k = 17). 

Other useful quantities to identify trapped modes are 
the rates of period changes (n/II) and the linear stabil- 
ity coefficients (77). Modes trapped in the envelope of the 
model should "feel" more strongly the effects of the surface 
gravitational contraction than untrapped modes, and thus 
the former should be characterized by lower values of H/H. 
This is clearly demonstrated in panel c of the figure, where 
we can see that the trapped modes are characterized by 
local minima in the distribution. On the other hand, it is 
well known from non-adiabatic arguments that the linear 
stability coefficients are larger for modes characterized by 
lower kinetic energies. This is depicted in panel d of the fig- 
ure, where the trapped modes (characterized by low kinetic 
energies) have local maxima in the redistribution. 

In view of the above discussion, since the mode at 
II m 468 s (k = 18) — which is identified as a trapped 
mode by FUEA07 — has a minimum in the observed and 
computed i?n-distributions, but it has a maximum in the 
kinetic energy, we conclude that this mode is not a trapped 
mode in the outer envelope. The mode at II w 401 s 
(k = 15), on the other hand, corresponds to a minimum in 
the observed and theoretical i?n-distributions, and a min- 
imum of the kinetic energy; so, we conclude that this is a 
genuine trapped mode in the envelope, confirming the con- 
clusion of FUEA07. However, the trapping cycle of about 
68 s (Afc « 3) suggested by FUEA07 is unvalidated in the 
frame of the present analysis since the mode at n m 468 s 
which is used by those authors would not be a trapped 
mode. 

In closing, a final note on the mode-trapping properties 
of our best-fit model is worth adding. The variations seen 
in the period distribution — as revealed by the i?n — n dia- 
gram of Fig. 5 — are due to mode-trapping effects inflicted 
by two chemical transition regions: the inner interface of 
O/C and the more external interface of O/C/He. The in- 
ternal chemical profile and the run of the Ledoux term B 
and the logarithm of the square of the Brunt- Vaisala fre- 
quency (N) of our best-fit model in terms of the outer frac- 
tional mass are depicted in Fig. 6. We can wonder at this 
point whether the O/C/He interface or the O/C chemi- 
cal transition region is more relevant at fixing the mode- 
trapping pattern of our model, or if there exists a sort of 
core-envelope degeneracy in the sense that both interfaces 
are equally effective in producing mode-trapping structure 
(see Montgomery et al. 2003). To gain some insight into 
this direction, we have redone our pulsation computations 
by minimizing the influence of a given chemical interface 



— amply used in asteroseismology — because of the many miss- 
ing modes in the period spectrum of PG 0122+200. 



8 The kinetic energy values correspond to a normalization of 
the radial eigenfunction of (, r /r = 1 at the stellar surface. 
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Fig. 6. The internal chemical profile of our best-fit model (T e g = 
81,540 K, Af, = 0.556M©) in terms of the outer fractional 
mass. Also shown are the profiles of the Ledoux term B and 
the logarithm of the square of the Brunt- Vaisala, frequency 
(N). The thickness of the outer envelope is M env = 0.024M Q . 
The location of the O/C and O/C/He chemical transition re- 
gions is emphasized with gray regions. The vertical line at 
— log(l — M r /M*) ~ 1.35 marks the location of the bottom 
of the envelope (Xh c ~ 0.14) [Color figure only available in the 
electronic version of the article]. 



on the period structure of the best-fit model 9 . We found 
that, at the domain of the observed range of periods in 
PG 0122+200, the mode-trapping features of our model are 
induced mostly by the chemical gradient at the O/C/He 
interface, being the O/C interface much less relevant. For 
periods longer than about 650 — 700 s, instead, it is the core 
chemical structure in the O/C interface that mostly fixes 
the mode trapping properties; this statement applies, for in- 
stance, to the cases of PG 1159-035 and RX J2117.1+3412 
(see Corsico & Althaus 2005, 2006 and Corsico et al. 2007 
for more details). 

4.5. The asteroseismological distance and parallax of 
PG 0122+200 

We employ the luminosity of our best-fit model to in- 
fer the seismic distance of PG 0122+200 from the Earth. 
First, we consider the flux predicted by a NLTE model 
atmosphere with T e g = 80 kK and logg = 7.5 inte- 
grated through the spectral response of the V filter. The 
model atmosphere was calculated with the Tubingen Model 
Atmosphere Package (see, for details, Werner et al. 2003). 
We obtain a bolometric correction BC= —5.89 and an 
absolute magnitude M v = 7.79. We account for the in- 
terstellar absorption, Ay, using the interstellar extinction 
model developed by Chen et al. (1998). We compute the 
seismic distance d according to the well-known relation: 
log d = g [m v — M v + 5 — Ay(d)] where the apparent mag- 
nitude is TO V = 16.8 ± 0.1 (FUEA07). The interstellar ab- 
sorption Ay(d) varies non linearly with the distance and 



We employ the same procedure like in Corsico & Althaus 
(2005, 2006); we refer the reader to those papers for details. 



also depends on the Galactic latitude (b) and longitude 
{I). For the equatorial coordinates of PG 0122+200 (Epoch 
B2000.00, a = l h 25 m 22 s .00, 6 = +20° 17' 54".0) the cor- 
responding Galactic coordinates are b = —41° 52' 1".2 and 
i = 133° 38' 16". 8. We solve for d and Ay iteratively and 
obtain a distance d — 614^2 P c an d an interstellar ex- 
tinction Ay = 0.0707. Note that our distance is w 13% 
smaller than the estimation of FUEA07 (d = 700±4°°°), 
and with its accuracy substantially improved. Finally, our 
calculations predict a parallax of tt ~ 1.6 mas. 

In closing, we estimate a "spectroscopic" distance of 
PG 0122+200. We first derive Ay = 3.1 E{B -V) = 0.19 
by employing the E(B — V) — 0.06 value from Dreizler 
& Heber (1998). The distance can be determined from the 
model y-flux comparing with my and using the spectro- 
scopic T e g and log g and the extinction Ay. We obtain a 
spectroscopic distance of 682 pc, and a parallax of ~ 1.47 
mas. We derive also an absolute magnitude M v = 7.55 
and a bolometric correction BC= —5.81 by employing the 
"spectroscopic" luminosity and radius — interpolated from 
the tracks assuming spectroscopic T e g and log g — and the 
flux predicted by the model atmosphere (see Table 2). 

5. Summary and conclusions 

In this paper we carried out an asteroseismological study 
of the cool pulsating PG1159 star PG 0122+200, a g-mode 
pulsator that defines the red edge of the GW Vir insta- 
bility domain at low luminosities. Our analysis is based 
on the full PG1159 evolutionary models of Althaus et al. 
(2005), Miller Bertolami & Althaus (2006) and Corsico et 
al. (2006). These models represent a solid basis to analyze 
the evolutionary and pulsational status of PG1159 stars 
like PG 0122+200. This is the second GW Vir star that is 
pulsationally analyzed in the frame of these state-of-the- 
art PG1159 evolutionary models — the first one being the 
hottest known GW Vir star, RX J2117. 1+3412; see Corsico 
et al. (2007). 

We first took advantage of the strong dependence of 
the period spacing of variable PG1159 stars on the stellar 
mass, and derived a value M t ~ O.625M by comparing 
All with the asymptotic period spacing of our models. 
We also compared All with the computed period spacing 
averaged over the period range observed in PG 0122+200, 
and derived a value of M* ~ 0.567M Q . Note that in both 
derivations of the stellar mass we made use of the spec- 
troscopic constraint that the effective temperature of the 
star should be ~ 80 kK. It is interesting to note that the 
stellar mass as inferred from the asymptotic period spacing 
is about 0.06M Q larger than that derived from the aver- 
age of the computed period spacings. This hints at possible 
systematics in the standard asteroseismological mass deter- 
minations methods, in particular when the full asymptotic 
regime (k 3> 1) has not been attained. We note that this 
systematics in the asteroseismological method is present 
not only in the case of full PG 1159 evolutionary models 
as we use here, but also in PG 1159 models artificially cre- 
ated (see Corsico & Althaus 2006). Because most analysis 
of pulsating PG1159 stars rely on the asymptotic period 
spacing, this point deserves to be explored for other GW 
Vir stars, issue which we address in a submitted paper. 

Next, we adopted a less conservative approach in which 
the individual observed pulsation periods alone — i.e., ig- 
noring "external" constraints such as the spectroscopic val- 
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ues of the surface gravity and effective temperature — natu- 
rally lead to an "asteroseismological" PG1159 model that is 
assumed to be representative of the target star. Specifically, 
the method consists in looking for the model that best re- 
produces the observed periods. The period fit was made on 
a grid of PG1159 models with a quite fine resolution in effec- 
tive temperature (AT e ff ^10 — 30 K) although admittedly 
coarse in stellar mass (AM* ~ 0.01 — 0.08M Q ). The match 
between the periods of the best-fit model and the observed 
periods in PG 0122+200 turns out be of an unprecedented 
quality for this type of studies, being the average of the pe- 
riod differences (observed versus theoretical) of only 0.88 s 
with a root-mean-square residual of 1.27 s. The stellar mass 
of the best-fit model is M* = 0.556±g;g?|M©. 

Interestingly enough, the mass of the best-fit model 
(M* = 0.556to'oi4-^©) i s closer to the spectroscopic value 
of M* = 0.53 ± 0.1M Q (Dreizler & Heber 1998; Miller 
Bertolami & Althaus 2006) than the asteroseismological 
mass derived in previous works, of 0.59 — 0.69M Q (FUEA07; 
O'Brien et al. 1998). 

Other characteristics of the best-fit model are summa- 
rized in Table 2. In particular, its effective temperature is 
nearly the same (to within 2%) as the spectroscopic T e ff. 
The surface gravity, on the other hand, is somewhat larger 
than the value given by spectroscopy. We also infer the 
"seismic distance" of PG 0122+200 by using the luminos- 
ity of our best-fit model. We obtain a distance d ~ 614 pc, 
somewhat smaller than that of FUEA07. 

Finally, our computations predict a temporal period 
drift for PG 0122+200 between 1.22 x 10~ 12 s/s and 3.26 x 
10~ 12 s/s. The positive values of II (increasing periods) re- 
flect the fact that our best-fit model is entering its white 
dwarf cooling domain where the effect of the increasing elec- 
tron degeneracy on the pulsation periods overwhelms that 
of the surface gravitational contraction, even for the modes 
trapped in the envelope. Strong theoretical arguments sug- 
gest that PG 0122+200 could be used to constrain the plas- 
mon neutrino rates in the dense interior of pre- white dwarfs 
on the basis of an observed value of II (O'Brien et al. 1998; 
O'Brien & Kawaler 2000). We defer a thorough exploration 
of this exciting issue to a forthcoming paper. 

The results of the period-fit procedure carried out in 
this work suggest that the asteroseismological mass of 
PG 0122+200 (~ 0.556M Q ) could be ~ 6 - 20% lower 
than thought hitherto (sec O'Brien et al. 1998 and more re- 
cently FUEA07) and in closer agreement (to within ~ 5%) 
with the spectroscopic mass derived by Miller Bertolami & 
Althaus (2006). This suggests that a reasonable consistency 
between the mass values obtained from both (very different) 
methods should be expected when detailed period-fit pro- 
cedures on full PG1159 evolutionary models such as those 
employed in this paper are considered. Even more, a bet- 
ter agreement between asteroseismological and spectroscopic 
masses of GW Vir stars could be found when the same evo- 
lutionary tracks are used for both the asteroseismological 
and the spectroscopic derivations of the stellar mass, as we 
do in the present work 10 . An anomalous case in this context 
could be RX J2117. 1+3412, for which we found an astero- 
seismological mass about 25% lower than the spectroscopic 
value by employing the same stellar evolution modeling 
than here (Corsico et al. 2007). As we suggested in that 



10 Sec Quirion et al. (2007) for an enlightening note about this 
topic. 



paper, the discrepancy in mass could be due to large errors 
in the spectroscopic determination of log g and T e ff for RX 
J2117. 1+3412, and/or uncertainties in the location of the 
evolutionary tracks in the HR and log T c ff — log g diagrams 
due to the modeling of PG1159 stars and their precursors 11 . 

In any case, detailed asteroseismological period-fits for 
other GW Vir stars based on full evolutionary models like 
we employ here, as well as precise spectroscopic determi- 
nation of the effective temperature and gravity of PG1159 
stars will be needed in the future if we want to reduce the 
persisting discrepancies in the stellar mass of these fasci- 
nating stars. 

In closing, in this paper we have been able to find 
a PG1159 model that nicely reproduces the period spec- 
trum observed in PG 0122+200 without artificially tuning 
the value of structural parameters such as the thickness 
of the outer envelope, the surface chemical abundances, or 
the shape of the core chemical profile which, instead, are 
kept fixed at the values predicted by our evolutionary com- 
putations. In some sense, this makes the fit derived more 
statistically significant. In particular, our PG 1159 evolu- 
tionary models are characterized by thick helium-rich en- 
velopes. However, we cannot discard the possibility that 
pulsating PG1159 star could harbor thin helium-rich en- 
velopes, a possibility sustained by the fact that PG 1159 
and born-again stars are observed to suffer from apprecia- 
ble mass loss. Resulting asterosesimological fits in this case 
would be worth exploring. 
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